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1 Introduction 

Determination of different biological factors on the base of dried blood spot 
technology has a great practical importance in investigation of back lands pop- 
ulations, in epidemiological studies or in special people contingents monitoring, 
see pQ, [2]. This technology presumes that blood sampling is performed by pa- 
tient himself, the sample is spotted on a dry, as a rule, porous surface (filter 
paper, cellulose acetate membrane, etc.), and the posterior transportation to a 
laboratory, for example, by post in a standard or a special envelope. 

Modern diagnostic equipment gives an opportunity to investigate many char- 
acteristics of dried blood spot, such as metabolites (see [3], [4], [5], [6], and [7]), 
hormones (see [8], [9], [10]), glycated hemoglobin (see [H]), and even some 
immune system parameters (see [H]). The possibility of DNA and RNA in- 
vestigations in such samples is of great importance, since it gives, for example, 
an opportunity for mass investigations of socially important infections such as 
AIDS, hepatitis, etc. 

In 1992, in USA a laboratory standard for dried blood spot testing (DBS) 
was elaborated, see [T3]. In 2001, in Russia, Application Instruction on Alcor Bio 
Ltd Reagent Kit for immuno-enzymatic determination of thyrotropic hormone 
in dried blood spot of newborn was approved. 

Under the dried blood spot technology using, the problem of the sample 
spotted volume determination remains one of the main practical questions. The 
existing versions of the technology may assume spotting of a known blood vol- 
ume by means of some dosing device and a posterior elution, or using of special 
filtering device for the plasma separation under dried spot preparation, see [14]. 

But there arc no any universal method calculating the volume of the blood 
spot, which does not use a dosing device. The solution of this problem gives an 
opportunity to increase essentially the accuracy of the results and to simplify the 
blood sampling procedure. There is a series of articles, where the calculations 
are based on the concentration of some electrolyte and on a correction of the 
plasma volume by the hematocrit value, see )16] . In the present work we try to 
elaborate a universal technology for the spotted volume calculations. 
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The solution of this problem is also important for such branches of Medicine 
as Catastrophe Medicine and Forensic Laboratory, where non-standard situa- 
tions are typical, and, for example, there is no opportunity to sample patient's 
blood, and so it is necessary to use the remains of the patient's blood on some 
other objects, instead of the standard blood sample. Such a situation can appear 
under investigations of the victims to traffic accidents, or other catastrophes. 

It is well-known that distinct biological indices (analytes) have distinct vari- 
ability, see [T5]. We try to use some mathematical algorithms to pick out a set 
of blood parameters which give an opportunity to retrieve the initial volume 
of the blood spotted, and use it to calculate exact concentrations of analyts 
interesting to a physician. For our analysis we used the database of biochemical 
blood parameters obtained in Russian Scientific Center of Roentgen-Radiology 
during 1995-2000, which includes more than 30000 of patients. 

2 Mathematical Model 

Let us describe the mathematical model of the problem. Let Xi, 1 < i < m, 
stand for the value of the result of the laboratory analysis on the ith molecular 
compound content. The value Xi obtained as a result of the blood sample 
analysis depends on two following parameters at least: the patient p which is 
selected from some collection ^f, and the volume A of the blood sample under 
the analysis. Thus, the value function of two parameters: 

Xi = Xi{p, A). 

The problem is to find a function y 2 , ■ ■ ■ , y m ) whose value 
/ {xi (P, A), x 2 (p, A), . . . , x m (p, A)) 

is close to A from the statistical point of view. 

Notice that due to the uniform distribution of the molecules under consider- 
ation in the blood, a fc-multiple extension of the volume must lead to the same 
increasing of all the indices Xi. In other words, we have Xi(p,kX) = kxi{p,X). 
Therefore, if / approximates the blood volume, then the following relation must 
be valid: 

f(xi(p,k\),...,x m (p,k\)) = f(kxi(p, A), . . . , fc x m (p, A)) w 

~ kf(x 1 (p,X), . . . , x m (p, A)). 

This notice is a natural motivation to look for the function / in the class of 
positively homogeneous functions of degree 1, i.e., we assume that the equality 

f(kyi,ky 2 , . . .,ky m ) = kf(y 1 ,y 2 , . . . ,y m ) 

holds for each positive k. Such functions are uniquely defined by their values 
at the unit sphere S"" -1 defined by the equation yf + ■ ■ ■ + yf n = 1. By g we 
denote the restriction of the function / onto this unit sphere. 
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Polynomials form the simplest but rich class of functions. Let us look for g 
among the functions which are the restrictions of the polynomials onto S" 1 ^ 1 . 
Our statistical experiments show that it is enough to consider the polynomials of 
degree two vanishing at the origin. In other words, wc put p — \/ y\ + ■ ■ ■ + y? n , 
and look for g in the form 

P rrr*. P P 
so the function / is supposed to be in the class 

/ = E a * Vi + E a ijViVilP' 

i i<j 

Thus, our problem is to find the coefficients a, and ctjj such that the func- 
tion obtained meets our objectives as well as possible. To formulate the latter 
condition mathematically, let us write down the following objective function. 

As we have already mentioned above, the available database gives us a table 
of specific values Xi S = x,(p s , A s ). We look for the function / such that the total 
squared deviation from the values A s is as small as possible. In other words, we 
have to find the on and Qy minimizing the objective function 

h = ^2(f(x ls , . . .,x ms ) - A s ) 2 = 

s 

= E (E x *> + E a ^ , + 3S , T 2 A -) • 

s V i i<j v Is r x rns / 

Notice that h considered as a function on on and aij is a non-negative 
quadric. In general position such a quadric possesses a unique minimum which 
can be found as a solution of linear equations system, i.e., from the condition 
that the differential of h vanishes. 

3 Application to the specific database 

The above algorithm determining the volume of a sample for calculation of the 
individual values of an arbitrary analyt was examined on the database of labo- 
ratory indices. The best results were obtained, when we reconstruct the volume 
by means of the following analyts: TP, K, Na. The correlation coefficients for 
the repaired and true values were 0.95-0.97. The algorithm obtained gives an 
opportunity to choose distinct sets of the indices for the volume reconstruction, 
that makes the algorithm multipurpose, i.e. it can be used for analysis of any 
laboratory blood indices. 

The method considered was applied to the specific database in RSCRR. 
This database was constructed from 35000 medical reports containing biochem- 
ical measuring data. We selected the reports containing the largest number of 
the biochemical data. So, we selected the set of 2637 cases with the next 17 
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biochemical data measured: Choi, TBil, DBil, TP, Alb, Urea, Crea, ALT, AST, 
Amy, ALP, K, Ca, Na, Fe, Glu, LDH. 

After calculation of the coefficients a, and o^- for the function /, we find 
out the following result: the number of patients p s which the inequality 

\f( x ls, ■ ■ ■ , %jns) — X s 

holds for, does not exceed 5%. This estimate agrees with the statistical signifi- 
cance of the result. 
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